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SMASHING THE GUITAR: AN EVOLVING NEUTRON STAR BOW SHOCK 
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ABSTRACT 

The Guitar nebula is a spectacular example of an Ha bow shock nebula produced by the interaction 
of a neutron star with its environment. The radio pulsar B2224+65 is traveling at ^800-1600 km s _1 
(for a distance of 1-2 kpc), placing it on the high- velocity tail of the pulsar velocity distribution. Here 
we report time evolution in the shape of the Guitar nebula, the first such observations for a bow shock 
nebula, as seen in Ha imaging with the Hubble Space Telescope. The morphology of the nebula provides 
no evidence for anisotropy in the pulsar wind, nor for fluctuations in the pulsar wind luminosity. The 
nebula shows morphological changes over two epochs spaced by seven years that imply the existence of 
significant gradients and inhomogeneities in the ambient interstellar medium. These observations offer 
astrophysically unique, in situ probes of length scales between 5 x 10~ 4 pc and 0.012 pc. Model fitting 
suggests that the nebula axis — and thus the three-dimensional velocity vector — lies within 20°of the 
plane of the sky, and also jointly constrains the distance to the neutron star and the ambient density. 

Subject headings: ISM:structure — pulsars:individual (PSR B2224+65) — shock waves — starsmeutron 
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1. INTRODUCTION 

While the steady decay of spin rates observed in radio 
pulsars provides a good estimate of the rate of energy loss 
E, only a small fraction of the energy output of a neutron 
star is typically converted to directly detectable electro- 
magnetic radiation. Most of the spindown energy loss is 
carried away by a relativistic wind, the properties of which 
are largely unknown. The best constraints on the relativis- 
tic wind derive from its interaction with the interstellar 
medium (ISM). This interaction has been observed, for 
example, in synchrotron nebulae such as the Crab Nebula 
(Kennel & Coroniti 1984; Gallant & Arons 1994; Melatos 
& Melrose 1996), where wisp structures are observed mov- 
ing at ~ 0.5c (Hester et al. 2002), as well as bow shocks 
observed in Ha emission from shock-excited neutral gas 
around PSR B1957+20 (Kulkarni & Hester 1988) and var- 
ious other pulsars (Chatterjee & Cordes 2002; Gaensler, 
Jones & Stappers 2002). 

The Guitar nebula was discovered in deep Ha imag- 
ing observations with the 5-m Hale telescope at Palomar 
Observatory (Cordes, Romani & Lundgren 1993). Us- 
ing a model for the Galactic electron density distribu- 
tion (Cordes & Lazio 2002), the dispersion measure (DM 

= f°n e ds = 35.3 pc cm~ 3 ) of PSR B2224+65 implies 
a distance of 1.9 kpc (hereafter parameterized as Di g). 
While the pulsar has a modest E= 10 33 1 erg s _1 , the large 
space velocity of the pulsar (~ 1640 Di g km s _1 ) provides 
the ram pressure needed to create a detectable bow shock 
nebula. The existence of an observable bow shock nebula 
also implies the presence of a significant neutral hydro- 
gen component in the ISM near the pulsar. Here we de- 
scribe Hubble Space Telescope (HST) observations of the 
time evolution of the nebula, and discuss the implications 
of the observations for our understanding of neutron star 
(NS) relativistic winds and the ISM. 

2. HST OBSERVATIONS AND MODELING 

High resolution HST observations obtained in 1994 De- 
cember with the Wide Field and Planetary Camera 2 



(WFPC2; Holtzman ct al. 1995) have been described pre- 
viously (Chatterjee & Cordes 2002). New WFPC2 obser- 
vations were obtained in 2001 December, with about 2.4 
times the exposure time of the earlier epoch. At both 
epochs, exposures were combined using variable-pixel lin- 
ear reconstruction ("Drizzling"; Fruchter & Hook 2002), 
yielding an effective pixel scale 0'.'0226, roughly half the 
size of the 0'.'0455 pixels on the WFPC2 Planetary Cam- 
era chip. The Ha images at the two epochs (Figure 1) 
were aligned to sub-pixel accuracy using eight stars (Fig- 
ure 2a). The tip of the nebula has moved 1.3" in 7 years, 
consistent with the radio proper motion of PSR B2224+65 
(fi = 182 ± 3 mas yr , at a position angle 52.1° ± 0.9°; 
Harrison, Lyne & Anderson 1993). However, unlike the 
simple translation with constant shape expected for a bow 
shock in a uniform ambient medium, the head of the neb- 
ula has changed morphologically, showing corrugations in 
the limb-brightened edge as well as gaps in the Ha flux, 
both of which vary with time. 

At the tip of the nebula, the stand-off radius of the bow 
shock has increased over the seven-year period, as deter- 
mined by model fitting (discussed below) at each epoch. 
Pressure balance between the relativistic NS wind and ram 
pressure from the ambient medium occurs at the stand-off 
radius, 

R = {E/^cnAmvD 1 / 2 , (1) 

where ua is the number density of particles with mean 
mass m and is the pulsar velocity. For the Guitar neb- 
ula, incorporating measured values for the pulsar E and 

i 1/2 

proper motion, the stand-off radius is Rq « 18 D 2 9 n A 
astronomical units (AU), and the predicted stand-off angle 
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milliarcseconds. Bright patches ap- 



pear along the edges of the nebula where it is constricted, 
as marked with arrows in Figure 1: the location of this 
brightening has moved ~ 0.8" parallel to the pulsar motion 
along both sides of the nebula. However, in Figure 2(a), 
the edges perpendicular to the nebula axis have moved 
outwards by less than 0.3", remaining essentially static at 
several points. Meanwhile, the rear edge of Ha emission 
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Fig. 1. — Ha images of the head of the Guitar nebula. The bottom panel shows a wide-field image of the Guitar nebula obtained at the 5-m 
Hale Telescope at Palomar (Chatterjee & Cordes 2002). High resolution images of the region marked with a box (~ 16" in size) were obtained 
with the HST PC in 1994 (right) and 2001 (left). North is upward and cast is to the left. Stars are circled, and arrows mark constricted 
regions in the limb-brightened H« emission that are conspicuously brighter (see text). Note that the neutron star is located at the very tip 
of the nebula at each epoch. 



has actually moved backwards by ~ 0.7" in 2001 com- 
pared to 1994. Possible explanations include turbulence 
in the shocked layers, their interaction with a complex 
magnetic field, or changes in any of the variables in the ex- 
pression for Rq (Equation 1), including E. However, even 
for the largest observed glitches (Hobbs et al. 2002), E 
changes temporarily by < 3%, a negligible variation com- 
pared to the large morphological changes observed. Be- 
sides being contrary to our general understanding of spin 
down, a time-variable E also fails to explain the bright- 
ening of the nebula at the locations where it appears con- 
stricted. Indeed, it is difficult to explain the evolving mor- 
phology without invoking variations in the ambient inter- 
stellar density, although instabilities in the shock structure 
are also possible (e.g. Draine & McKee 1993). We propose 
that the observations reflect the motion of PSR B2224+65 
through random density inhomogeneities combined with a 
gradient towards a region of lower density. The bright rear 
edge (~ 11" downstream of the nebula tip in 1994) marks 
a sharp increase in density which the NS broke through 
~ 70 years ago. Presently, the shock at the rear confines 



the relativistic NS wind, leading to a brightening of the 
head of the nebula and preventing the wind from power- 
ing the rest of the Guitar body, thus producing the dim, 
elongated and narrow "neck" of the Guitar in Figure 1. 

We argue further that the morphology of the head of 
the Guitar nebula is an analog for confinement by another 
high-density region ~ 300 years ago, which created the 
rounded end of the Guitar body. Fluctuations in the am- 
bient density cause the body to be brighter where it ap- 
pears constricted. In a few hundred years, as the larger 
Guitar body fades from view, what is currently the head 
of the nebula may expand to become another guitar-like 
structure. The scenario described here, while plausible, 
needs to be verified with future high-resolution monitor- 
ing observations as well as time-dependent hydrodynamic 
modeling of shock fronts in an ambient medium with sig- 
nificant density fluctuations. 

Currently, in order to quantify the change in stand-off 
radius between epochs, we have modeled the shock front 
under the assumption that the nebula tip is in quasi-static 
equilibrium with the ambient medium. A momentum- 
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Table 1 

Best-fit Model Parameters for HST Images 





Position Angle 


Inclination''' 


Scale Factor 


e 


Epoch 


(°) 


(°) 


p k - p W /2 ) 


(") 


1994 


48 ±2 


90 ±30 


3.6 ±1.2 


0.12 ±0.04 


2001 


50 ±2 


90 ±20 


4.4 ±1.2 


0.15 ±0.04 



^The fit for inclination angle was inconclusive, especially for the 
1994 data, as discussed in Chatterjee & Cordes (2002). 
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Fig. 2. — (a) Left: The outline of the limb-brightened head of 
the Guitar nebula in 1994 (filled circles) and 2001 (open circles), 
obtained by eye. The images are aligned using 8 bright stars in 
the PC image at each epoch: to demonstrate the image registra- 
tion, two stars in the PC image section displayed in Figure 1 are 
marked with crosses in 1994 and open circles in 2001. (b) Right: 
Slices through the best fit point on the \ 2 surface for the model fit 
procedure on 1994 data (filled circles) and 2001 data (open circles), 
showing the relative change in x 2 as a function of the scale factor 

'kpc n A 

nebula tip is 71,4 (2001)/nA (1994) 0.7. 



conserving bow shock model (Wilkin 1996) was adapted 
to fit the Ha emission at each epoch. The momentum- 
conserving description has known limitations (Bucciantini 
& Bandiera 2001), especially since it applies only to an 
ambient medium of uniform density. To avoid these prob- 
lems, we restricted the model fit to within 2'.'6 of the tip of 
the nebula, where it is smooth and symmetric. The model 
is parameterized by the position angle and the inclination 
of the nebula to the line of sight (neither of which are ex- 
pected to change significantly between epochs), the thick- 
ness of the shocked layer that emits Ha, and by a scale 



factor S = D, 



-1/2 



y kpc' L A > which isolates the dependence of 
the apparent size of the nebula on distance and ambient 
density. Details of the model fitting procedure (for the 
1994 data) are given in Chatterjee & Cordes (2002) and 
the best-fit parameters for both epochs (position angle, 
inclination, scale factor and stand-off angle) are listed in 
Table 1. At both epochs, the fit constrains the nebula to 
lie in the plane of the sky. As shown in Figure 2(b), the 
best-fit scale factor varies significantly between 1994 and 
2001. Since the fractional change in distance to the nebula 
is negligible over seven years, the change in scale factor 
implies a decrease in ambient density by a factor ~ 0.7 
(from 0.006 cm~ 3 to 0.004 cm~ 3 for D = 1.9 kpc) over 
1.3", corresponding to a length scale of 2500 Di. 9 AU. The 
implied change in DM for density changes on this length 
scale is ~ 10~ 4 pc cm~ 3 , which may be detectable with 
sensitive pulse timing observations. Additionally, the fits 



establish a joint constraint on the distance to the NS and 

the ambient density, Dk pc = OASn^^ 4 . We note that 
the densities obtained above are low, suggesting a possible 
overestimate of the distance. For an ambient density of 
0.05 cm~ 3 , which is comparable to the density of the lo- 
cal warm ionized medium (e.g. Paresce 1984), the implied 
distance is 1 kpc (and the height above the Galactic plane 
is reduced from 240 pc to 120 pc). 

To check for unresolved small-scale structure in the 
limb-brightened nebula, the autocorrelation function was 
calculated for sections of the image with and without neb- 
ular emission. After accounting for Poisson noise and the 
contribution from smooth extended structure, the excess 
in the on-nebula autocorrelation function due to barely- 
resolved or unresolved structure (50 mas or less) is < 5% 
of the larger-scale nebular emission. We conclude that 
50 mas represents a lower limit on the angular scale of 
structure in the interstellar medium probed by these HST 
observations of the Guitar nebula. 

3. DISCUSSION 

The wavenumber spectrum for electron density fluctu- 
ations in the general ISM has been delineated through 
a variety of measurements, including radio scintillation, 
scattering, pulse time-of-arrival, and Faraday-rotation, as 
illustrated in Figure 3. These measurements constrain the 
power levels in different wavenumber intervals, and sug- 
gest an overall consistency with a Kolmogorov turbulence 
process (Armstrong, Rickett, & Spangler 1995), although 
this consistency is only coarse and may be an illusion as- 
sociated with the large amplitude scale of the diagram. 
Departures from the Kolmogorov spectrum are also in- 
ferred from pulsar scintillation measurements. The HST 
measurements of the Guitar Nebula provide constraints 
on wavenumbers not easily accessible by other probes and 
therefore provide a new tool for investigating fine structure 
in the interstellar density. 

Along with the properties of the ISM, the Guitar neb- 
ula also provides a probe of the pulsar's relativistic wind. 
The properties of pulsar winds, including the magnetiza- 
tion parameter a, the ratio of the Poynting flux to the 
kinetic energy flux, have been inferred primarily from the 
Crab nebula (Kennel & Coroniti 1984; Gallant & Arons 
1994; Melatos & Melrose 1996). The existence of bow 
shock nebulae such as the Guitar requires collisional exci- 
tation of the neutral interstellar medium through interac- 
tions with shocked electrons and protons. These charged 
particles can be ejected from the NS itself, or originate 
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Fig. 3. — Power spectrum for electron density variations in the 
ISM, defined such that the integral over all wavenumbers is the 
mean-square electron density. The region labeled 'GN' designates 
the constraints derived from HST observations of the Guitar Nebula 
described here. The wavenumber range corresponds to the recipro- 
cals of the length scales probed on angular scales between ~ 0.05 
and 1.3 arcsec, assuming a distance of 1.9 kpc and a combined filling 
factor and ionized fraction ranging from 0.02 to 0.5. The dashed line 
labelled 'HI' is an estimate using the length scales of standard HI 
clouds (McKee & Ostriker 1977), assumed to be partially ionized. 
The line labelled 'MSP' represents variations in dispersion measure 
for the millisecond pulsar B1937+21 (Cordes ct al. 1990; Kaspi, Tay- 
lor, and Ryba 1994) The observational constraints at the highest 
wavenumbers result from measurements of interstellar scintillation 
and scattering. The point for Sgr A* represents the power level for 
scattering of the Galactic center source, corrected for the proxim- 
ity of the scattering screen to the source (Lazio & Cordes 1998). 
The hatched regions incorporate constraints from scintillation and 
scattering observed for several hundred sources (NE2001; Cordes & 
Lazio 2002, and references therein). The light solid line designates 
a power spectrum with power-law index equal to the Kolmogorov 
value (-11/3). 



from interstellar atoms through photoionization or mag- 
netic reconncction outside the pulsar light cylinder radius 
(R LC = cP/2ir). The process by which the NS Poynt- 
ing flux is converted to particle kinetic energy is not well 
understood: overviews of different processes are provided 
by Begelman & Li (1992), Gallant et al. (2002, ion load- 
ing) and Lyubarsky & Kirk (2001, striped pulsar winds). 
At the stand-off radius (Rq ~> Rlc), pressure balance re- 
quires that the spindown energy of the NS be carried by 
the particle flux (a -C 1), but the shape of the nebula may 
encode information about the (rotation-averaged) shape of 
the NS wind. 

Along with the Guitar nebula, future observations of 
the evolution of bow shock nebulae will be possible for the 
nearby radio-quiet neutron star RX J1856.5— 3754 (van 
Kcrkwijk & Kulkarni 2001), only ^120 pc away (Ka- 
plan, van Kerkwijk & Anderson 2002; Walter & Lattimer 
2002). Even at its relatively low speed (< 200 km s _1 ), 
the neutron star travels ~ 1" in a year, and evolution of 
the nebula should be evident on roughly this time scale. 



PSR J2124— 3358, a nearby millisecond pulsar with a com- 
plex bow shock nebula (Gacnslcr, Jones & Stappcrs 2002), 
is also promising in this regard, while the discovery of a 
bow shock nebula powered by the Poynting flux from a 
magnetar would allow investigation of the relativistic wind 
in the presence of an ultra-strong magnetic field. 

This work is based in part on observations made with 
the NASA/ESA Hubble Space Telescope, obtained at the 
Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observa- 
tions are associated with proposals 5387 and 9129. This 
work was also supported by NSF grants AST 9819931 
and AST 0206036, and made use of NASA's Astrophysics 
Data System Abstract Service and the arXiv . org astro-ph 
preprint service. 
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